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Abstract—In this work, we employ a novel basis matrix 
method, applied in conjunction with the singular value 
decomposition (SVD) algorithm to achieve sub-wavelength 
resolution in the images derived by using phase conjugation 
(PC) scheme in the microwave imaging problem. [luminative 
examples of different bar-code type of distributions are 
presented, together with the images recovered by using the 
proposed method. 
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I INTRODUCTION 


Recently, the topics of phase conjugation (PC) [1] as 
well as time-reversal mirrors (TRM) [2] have been 
investigated in the context of microwave imaging. A 
number of papers that deal with the methods based on far- 
field measurements have been reported in the literature, 
including: those that employ metallic strip gratings which 
perform evanescent-to-propagating wave conversion [3]; 
and split-ring resonators loaded with varactor diodes [4]. 
The resolution capability of a phase conjugating lens has 
also been studied in [5]-[6]. In this work, we propose a basis 
matrix method applied in conjunction with the singular 
value decomposition (SVD) to achieve a sub-wavelength 
level of resolution in the context of microwave imaging 
problem using phase conjugation lens. We described the 
proposed method in detail and examined its abilities to 
resolve at the sub-wavelength level. 


II. IMAGING ALGORITHMS 


A. Phase Conjugation 


Figure 1 shows a representative scheme for 3D phase 
conjugation lens imaging system. The test objects are 
perfect electric conductor (PEC) squares that are distributed 
in the source plane (f(x,y)@z=0), and are illuminated by an 
incident plane wave. At the measurement plane 
(g(x,y)@z=d), the phase conjugation lens collects the 
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scattered field, phase conjugates this field and propagates it 
to the image plane (A(x,y)@z=2d) to form the image. 


For the purpose of our simulation, we mimic the phase 
conjugating operation of the lens, applied to the measured 
field, on the computer. Let us assume that we are working 
with the £,-component of the electric field. Then, the steps 
involved in the phase conjugation procedure can be 
summarized as follows: (i) measure the £,-component of the 
scattered electric field from the PEC squares at the 
measurement plane, which is located at z=d; (ii) phase 
conjugate and Fourier transform the spatial distribution to 
the wavenumber domain; (iii) propagate the spectral 
distribution to the image plane at z=2d; (iv) inverse-Fourier 
transform the spectral distribution at z=2d to obtain the 
image. 
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Fig. 1. Schematic for 3D imaging via basis matrix method. 


Note that in the proposed method we only need to deal 
with the data sampled in the image plane to improve the 
resolution of the image. This, in turn, results in considerable 
saving in terms of data collection and processing, in 
comparison to that needed if we’re to collect the data 
directly in the measurement plane. Note also that the phase 
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conjugation can be done by using a lens placed at the 
measurement plane and this lens can be flat, if desired. 


B. Basis Matrix Method with Singular Value 
Decomposition 


Next, we describe an imaging algorithm based on the 
basis matrix method. As a starter, we use a MoM code 
(FEKO) to simulate the case of a single PEC square at a 
time, which is located in the object plane. We derive the 
scattered field distribution in the measurement plane, phase 
conjugate it, and then let it propagate to the image plane. 
Finally, we extract a part of the image distribution, for 
instance its values along the y=0 cut to form a column of the 
basis matrix. The number of the columns in this matrix 
corresponds to the number of possible positions for the PEC 
square in the source plane. We apply the singular value 
decomposition algorithm to the basis matrix, and use a 
threshold to delete the singular vectors, whose singular 
value fall below the threshold. 


Ill. NUMERICAL RESULTS 


Figure | shows a representative scheme 3D imaging by 
using the proposed method. The size of the PEC square is 
0.052*0.052, and the spacing is 0.252 in both x- and y- 
directions. The distance between the source and 
measurement planes is 34. The size of the measurement 
aperture is 301*30A with 0.054 interval in both x- and y- 
directions. And the size of the image zone is /04*/0A with 
0.254 interval in both x- and y-directions. The operating 
frequency is /GHz, and the plane wave is normally incident. 
Fig. 2 shows the imaging models for the 1D and 2D cases. 
The basis matrices are comprised of /3 and 39 columns, 
respectively. The object zone in the x-direction ranges from 
-1.54 to 1.52. 
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Fig. 2. Imaging models for (a) 1D case (x-direction); and (b) 
2D case (x-y plane). 
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Figures 3 and 4 show the recovered results for the 1D 
and 2D cases, respectively, the stars and the bars represent 
the original and the recovered positions of PEC squares, 
respectively, and the horizontal line at the level of 0.5 
(normalized) is used to distinguish the real objects from the 
false ones. Here, we have set a threshold of /0° for the 
truncation of the singular values. For the 1D case, we only 
extract y=0 cut to generate the basis matrix, and find /3 
singular vectors survive. For the 2D case, we use three cuts 
e.g., at y=-0.1A, OA, 0.12, and 35 singular vectors are 
survived. From Figs. 3 and 4, we see that the basis matrix 
method, applied in conjunction with the SVD is capable of 
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accurately recovering the position of the 


systematically investigated the noise. 


object 
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Fig. 3. Imaging results for the 1D case (x-direction) via the 
basis matrix method. 
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Fig. 4. Imaging results for the 2D case (x-y plane) via the 
basis matrix method. 


We have considered the noise issue, and have also 
studied the 2D case in the x-z plane, i.e., objects with depth 
information. These results will be included in the 
presentation. 
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